Structural signature of jamming transition in thermal amorphous systems 
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In thermal amorphous systems, the first peak of the pair correlation function g(r) shows a maxi- 
mum height (?i lax at a volume fraction cj> = cj> v that increases with the temperature, gj 13 * diverges 
at the T = jamming transition at 4> = 4> c . Molecular dynamics simulations show that some typi- 
cal quantities, such as the pressure, bulk modulus, shear modulus, and boson peak frequency that 
behave power law scaiings with — C in marginally jammed solids at T — 0, all show scaiings with 
4> — 4>c when tf> > <f> v , while the scaiings break down when <f> < cj> v . The presence of g^ ax is thus 
not only a thermal vestige of the T — jamming transition, but more importantly the structural 
signature of the jamming transition. 

PACS numbers: 63.50.Lm, 64.70.Pv, 64.70.Pm 



Amorphous systems undergo jamming (liquid-like to 
solid-like) transition via lowering the temperature T, in- 
creasing the volume fraction 0, or decreasing the shear 
stress <J xy , as described by Fig. [T]of the jamming phase 
diagram [l|, 0, At T — and a xy = 0, the jamming 
transition occurs at a critical volume fraction <j) c , labeled 
as Point J in Fig. [TJ which is the jamming transition 
of packings of frictionless spheres interacting via repul- 
sions. Point J shows criticaiitv with multiple scaiings and 
diverging length scales J, i, 0, S 1 M M M M M ■ 
Being the only singular point in the jamming phase di- 
agram, Point J may control the jamming transition at 
T > and a xy ^ and the properties of jammed solids. 
As shown in recent studies on boson peak in realistic 
glasses 0,[I1] and energy transport in jammed solids [3], 
properties of low temperature glasses can be understood 
from marginally jammed solids near Point J. 

Approaching the jamming transition from the un- 
jammed side in the T — <fi plane of the jamming phase 
diagram, the viscosity or the relaxation time of the liquid 
grows up rapidly and exceeds the measurable time win- 
dow at the glass transition. It is thus natural to adopt 
the glass transition as the jamming transition at T > 0. 
However, the glass transition is quantitatively arbitrary, 
because ideally the glass transition corresponds to infinite 
relaxation time which is inaccessible and no structural 
signature of the glass transition has ever been observed. 
The ambiguity of the glass transition thus weakens its 
validity as the candidate to the jamming transition. 

It has been recently reported that, when the volume 
fraction of an amorphous system is varied at a fixed 
temperature, the height of the first peak of the pair 
correlation function g(r), g\ (illustrated in the inset to 
Fig. H[a)), reaches the maximum value gf 1 ^ at a tem- 
perature dependent volume fraction <j> v [lj. This struc- 
tural maximum is the thermal vestige of the T = jam- 
ming transition at Point J at which g™ ax — > oo [16] . The 
occurrence of g™ ax is selective about trajectories: g\ in- 
creases monotonically with decreasing the temperature if 
the volume fraction is fixed, so g , 5 nax exists only when the 
volume fraction varies. Computer simulations have also 
shown that the occurrence of g , i tlax is not consistent with 




FIG. 1: (Color online) Schematic plot of the jamming phase 
diagram determined by the temperature T, volume fraction <j>, 
and shear stress a xy . Systems in the blue region are jammed. 
Point J is the jamming transition at T = and o xy — 0. 



the glass transition determined by the relaxation time 
exceeding the measurable time window [1]. Therefore, 
<7™ ax is not a signature of the glass transition. 

In this letter, we measure material properties of ther- 
mal amorphous systems in the vicinity of <p v by tuning 
the volume fraction <f> at fixed temperature. The pressure 
P, bulk modulus B, shear modulus G, and boson peak 
frequency lo* , which show power law scaiings with 4> — 4> c 
in marginally jammed solids at T = [H, [y], are scaled 
well with <p — 4> c when <fi > cf> v . The scaiings break down 
when <p < <p v . The presence of g" lax is thus not merely 
the thermal vestige of the T = jamming transition. It is 
a robust structural signature of the jamming transition. 

We perform molecular dynamics simulations at con- 
stant NPT in a three-dimensional cubic box with peri- 
odic boundary conditions in all directions 0, [l7| ■ The 
length of the box L varies in time t to maintain con- 
stant pressure. A 50 : 50 binary mixture of N = 10, 000 
particles with equal mass m move in the box under the 
constraint of the thermostat and barostat [HI]. The di- 
ameters of the small and large particles are a and 1.4<r, 
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respectively. Particles i and j interact via a short range 
repulsive potential: V^r^) = e(l — / Oij) a / a when 
their separation is less than the sum of their radii 
<7jj = ((Tj + Oj)/2, and zero otherwise. We use a = 2 
(harmonic) and 5/2 (Hertzian) to study the potential de- 
penden ce. The length, time, and energy are in the units 
of <7, v/mcr 2 /e, and e. The Boltzmann constant ks is set 
to unity. 

We measure the pair correlation function g(r) for large- 
large particles. As shown in Fig. H{a-b), for both har- 
monic and Hertzian repulsions, g±, the height of the first 
peak of g(r), is maximal at = <p v . In the range of 
temperatures studied, Fig. [3ja) shows that 

T«(0„-0 c ) a , (i) 

where C » 0.650 ± 0.002 Eq. Q arises from the 

relaxation of the bond length between interacting 
particles under thermal perturbations As shown in 
Fig. EJb) , the peak value <jrj nax increases with decreasing 
the temperature and diverges at T = (Point J) accord- 
ing to the power law 

ff r x «(^-0c)~ - 9 , (2) 

which is independent on the particle interactions. The 
dependence of (?f lax on <f> v is the same as gi(4>) in 
marginally jammed solids at T = [l6l |. which hints that 
gm&x comes a l ng with a typical material property change 
presented in jammed solids. The value of g™ ax does not 
show any apparent system size or history dependence [l[ . 

At T = 0, multiple quantities such as the pressure, 
bulk modulus, shear modulus, boson peak frequency, and 
coordination number (number of interacting neighbors 
per particle) show power law scalings with <j> — C in 
marginally jammed solids d,@|. The pressure P is scaled 
with ((j)-(j) c ) a - 1 @. At T > 0, we calculate the pressure 
from the virial expression P = (NkpT + £V F^ )/L 3 , 
where Fij — —dV(rij)/drij is the force between parti- 
cles i and j and the sum is over all pairs of interacting 
particles. When the kinetic energy is much lower than 
the potential energy, fluctuates around its equilibrium 
T = value r^-. The fluctuations are so small that they 
only contribute a tiny correction to the T = pressure. 
We thus expect that P is still proportional to (0 — C ) Q ~ 1 
at T > if the systems are jammed. Fig. [U[c-d) show 
that when > 0„, 

p = A P (4>-<j )c r- 1 , (3) 

where Ap is nearly insensitive to the temperature; when 
4> < 4> v , Eq. stops working. Accordingly, when > 
4> v , there is a temperature independent bulk modulus 
B — on (0 — c ) a ~ 2 , while the bulk modulus varies 
with both volume fraction and temperature when < 4> v . 

At T = 0, the shear modulus G is scaled with 
(0 — (f> c ) a ~ 3 ^ 2 in marginally jammed solids 5]. The 
ratio of the moduli G/B oc (0 — 0c) 1 / 2 that decays 
to zero at Point J, while in normal solids like crys- 
tals, G/B does not show strong volume fraction depen- 
dence. At T > 0, we measure the storage shear modulus 




FIG. 2: (Color online) The height of the first peak of g{r), gi, 
pressure P, storage shear modulus G' , and boson peak fre- 
quency lo* , versus the volume fraction — C for systems of 
N = 10, 000 particles interacting via harmonic (left column) 
and Hertzian (right column) repulsions at T = 10~ J . The 
inset to (a) shows the pair correlation function g(r) measured 
when gi reaches its maximum. The red dashed line illus- 
trates the definition of gi. The red solid curves are (c) P = 
0.24(0 - C ), (d) P = 0.18(0 - C ) 3/2 , (e) G' = 1.04(0 - C ), 
(f) G' = 0.55(0 - 4> c f /2 , (g) uj* = 0.31(0 - C ) 1/2 , and (h) 
lj* — 0.22(0 — </> c ) 3,/4 , respectively. The frequencies at which 
G' were measured are 0.01 (blue squares), 0.001 (magenta 
triangles), and 0.0001 (black circles), respectively. The ver- 
tical blue dashed lines are guide to eye, which denote where 
gi reaches its maximum. The critical volume fraction of the 
T = jamming transition <f> c is 0.650 ± 0.002. 



G' at various frequencies to by applying an oscillatory 
shear strain jsm(ujt) in the x— direction (shear gradient 
is in the y— direction) and measuring the shear stress 

a xy = -(I2iPx,iPy^/ m + I2i3 X ij F y^j)/ L3 ' where p Xfi 
and p y> i are the x— and y— components of the momen- 
tum of particle i, Xij is the ie— component of r^, and 
F Vt ij is the y— component of We measure G' at suf- 
ficiently small 7, so that linear response is guaranteed. 
Again, when the kinetic energy is negligible compared to 
the potential energy, the shear stress is determined by 
the particle interactions. In Fig.[^e-f), we show G' mea- 
sured over two decades of low frequencies. In the range 
of frequencies and volume fractions studied, G' does not 



3 




FIG. 3: (Color online) (a) Temperature T and (b) maximum 
value of gi, gi lax 5 versus <j> v — <j)c for harmonic (black circles) 
and Hertzian (red squares) repulsions, where <j) c — 0.650 ± 
0.002. The black solid and red dashed lines have a slope of 2 
and 5/2 in (a), and —0.9 in (b). 

depend on the frequency. When <fi > (fi v , 

G' = A G (T){<t>-4> c ) a '\ (4) 

where Aq increases with decreasing the temperature. 
When cf> < cf> v , Eq. (j4|) no longer holds. The temper- 
ature dependence of A G is not surprising. At a fixed 
volume fraction <fi > <fi v , if we increase the tempera- 
ture of a jammed solid approaching unjamming transi- 
tion, the shear modulus G'{T) must decrease. A G — 
G'(T)/(cj) — 4> c ) a " 1 thus decreases with increasing the 
temperature. 

The scaling of G' with (j> — <fi c is different from that 
in marginally jammed solids at T = 0. This discrep- 
ancy should arise from the thermal motion of particles 
around their equilibrium positions. Unlike in crystals 
in which the shear modulus G is determined by the re- 
sponse of the inter-particle forces to the strain, e.g. G 
is almost constant upon compression for harmonic re- 
pulsions (a = 2), in jammed systems the heterogeneous 
inter-particle bond orientation and forces are coupled to- 
gether, which lead to a much more complicated volume 
fraction dependence of the shear modulus. The thermal 
motion of particles increases the heterogeneity, so that 
the volume fraction dependence of the shear modulus at 
T > is pushed further away from crystals. The ratio of 
the moduli G/B oc (0— <fi c ) shows a stronger volume frac- 
tion dependence than that in marginally jammed solids. 
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FIG. 4: (Color online) (a) Density of vibrational states D(lu) 
and (b) reduced density of vibrational states D{uj)/uj 2 of 
systems with harmonic repulsions at T = 10~ 5 and <f) — 
0.6843 (black solid), 0.6760 (black dashed), 0.6696 (black dot- 
dashed), 0.6606 (red dotted), 0.6535 (blue solid), and 0.6466 
(blue dashed). The red dotted curves are measured approx- 
imately at 4> v at which the height of the first peak of g(r) 
reaches the maximum. 

The vibrational properties of marginally jammed solids 
have shown some anomalies. For instance, in marginally 
jammed solids, there is a plateau in the density of vibra- 
tional states D{u>) at low frequencies [f|. The plateau 
extends to lower frequencies approaching the unjamming 
transition. Excess low-frequency vibrational modes be- 
yond Debye's prediction aggregate in the vicinity of a 
characteristic frequency uj* and build up the boson peak 
[20l |2lj . The boson peak frequency ui* is scaled with 
{(j)-^)^-^' 2 [SHGj. At T > 0, the density of vibra- 
tional states can be calculated from the Fourier transform 
of the velocity correlation C(t) = (v(t) ■ v(0)) (22J: 

D{u) = -^J C(t)cos(ujt)dt, (5) 

where (.) denotes the average over particles and ensem- 
bles. Fig. IDJa) shows the density of vibrational states of 
systems with harmonic repulsions measured at T = 10~ 5 
and various <f>. When <j) > cf> v , a plateau is developed in 
D(lu). The plateau shifts to lower frequencies and grows 
up with decreasing <f>. When <f> < <f> v , the plateau dis- 
appears and is replaced with a peak that captures the 
general feature of the density of states of liquids [22l |. 
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The disappearance of the plateau near <p v is a sign of the 
unjamming transition. 

In Fig.[Uh) we highlight the boson peak by plotting the 
reduced density of vibrational states D(uj)/lu 2 . The bo- 
son peak shifts to lower frequencies, while its magnitude 
increases, with decreasing (f>. As shown in Fig. [2jg-h), 
jj° ax signifies the crossover between two regimes of dif- 
ferent volume fraction dependence of w*. When <f> < <f> v , 

u*=A u {T){<i>-<l>c) {a - 1)/ \ (6) 

where A u decreases with increasing the temperature. 
Eq. ([6]) starts to break down when <j) < <p v . The tem- 
perature dependence of is understandable from the 
fact that the boson peak in glasses shifts to lower fre- 
quencies with increasing the temperature [13, HH ■ At a 
fixed 4> > 4>v, = lj* /(4> - </>c) < ' Q ~ 1 ^ 2 thus decreases 
with increasing the temperature. 

In marginally jammed solids, another important quan- 
tity that shows the power law scaling with 4> — <p c is 
the coordination number z: z — z c oc {<f> — 0c) 1 / 2 . At 
Point J, z = z c = 2d is the isostatic value, the minimum 
constraint requirement to maintain mechanical stability, 
where d is the dimension of space. In thermal systems, 
z fluctuates in time, which is no longer a valid criterion 
of mechanical stability. Although the time averaged co- 
ordination number drops quickly when cf> < <j) v , the coor- 
dination number does not show any well-defined scalings 
with <fi — 4> c . 

The scalings shown in Eqs. ©-(HI) and ([6j) demonstrate 
that thermal amorphous systems at 4> > <j) v are jammed, 
since similar scalings exist in marginally jammed solids 
at T = 0. The temperature dependent crossover volume 
fraction (f> v separates two regimes with apparently differ- 



ent material properties. The occurrence of the maximum 
height of the first peak of g(r) is thus a structural signa- 
ture of the jamming transition in general. For systems 
studied, the jamming transition in the T — (f> plane of the 
jamming phase diagram shown in Fig. |T| is well-defined 
by Eq. (p}, which converges to Point J at T = 0. 

Eqs. (P)-(|ll) and © also highlight the key role of Point 
J in controlling the jamming transition and material 
properties of jammed solids. Point J is the only critical 
point in the T — <p plane of the jamming phase diagram, 
at which diverging length scales are able to exist. During 
the jamming transition of thermal amorphous systems we 
do not expect to observe diverging length scales because 
the critical scalings are all controlled by <f> c . 

For thermal systems as shown in Fig. [2 the relaxation 
time exceeds the measurable time window and shows 
rapid divergence at a much lower volume fraction than 
<p v . The jamming transition signified by is thus 

deeply inside the glass transition at which the system is 
already a glass. It is ambiguous whether there exists the 
thermodynamic glass transition and whether the colloidal 
glass transition in the hard sphere limit happens at Point 
J. The jamming transition interpreted by the glass tran- 
sition is thus not well-defined. Our study propose an un- 
ambiguous liquid-solid transition of amorphous systems 
determined by a robust structural signature, which opens 
a new point of view to understand the phase transition 
of amorphous systems and inspires new experiments [23| . 
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